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We review the electromagnetic form factor of heavy quarks with empha- 
sis on the QCD radiative corrections at two-loop order in the perturbative 
expansion. We discuss important properties of the heavy-quark form fac- 
tor such as its exponentiation in the high-energy limit and its role in QCD 
factorization theorems for massive n-parton amplitudes. 
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1. Introduction 

The electromagnetic form factor of heavy quarks comprises the simplest 
example of a scattering amplitude in Quantum Chromodynamics (QCD) 
for the study of mass effects in QCD hard scattering processes including 
radiative corrections at higher orders. As such, it has received considerable 
attention in recent years. 

Explicit results for the two-loop QCD corrections have been obtained 
for the vector- and axial-vector coupling as well as the anomaly contribu- 
tions [1—3]. Moreover, universal features of higher order radiative correc- 
tions, such as the exponentiation of logarithms in the heavy-quark mass 
m together with the infrared (soft) singularities within dimensional reg- 
ularization have been addressed [4]. As a consequence, these factorization 
properties in the soft and collinear limit imply for on-shell massive n-parton 
QCD amplitudes a simple relation between massless and massive scattering 
amplitudes at any order in perturbation theory [4,5]. 

In the following, we review briefly recent progress in studies of QCD 
corrections to the heavy-quark form factor. 
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2. QCD radiative corrections 

The object of our interest are the scalar functions T\ and T 2 which 
parametrize the vector current of an on-shell heavy-quark pair of momenta 
k\ and k 2 and of mass m. We denote the associated vertex function as 
and for a photon of (space-like) virtuality Q 2 = —q 2 > 0, we can write 

r„(fci,fc2) = (i) 

ie q ^(fc 1 ) ^ J r i(Q 2 ,m 2 ,a s ) - -^cr^ u q u J r 2 (Q 2 ,m 2 ,a s )j ip(k 2 ) ■ 

Here e q is the charge and T\ and T 2 are the electric and magnetic space-like 
quark form factors. They are gauge invariant, but in perturbative QCD at 
higher orders in general divergent and their power expansion in the strong 
coupling a s at the scale fj? reads, 



where we have introduced the shorthand notation a s = a s /(47r). 

Up to the two-loop level there is only a relatively small number of Feyn- 
man diagrams to be computed. All divergences are regulated in D = 4 — 2e 
dimensions and the results are given as a Laurent expansion in e. Except 
for additional diagrams from heavy-quark self-energy corrections on the ex- 
ternal lines, which contribute to the wave function renormalization, the 
diagrams are in one-to-one correspondence to the calculation of the form 
factor of massless quarks (see e.g. Ref. [6]). 

The reduction of all Feynman integrals to so-called master integrals pro- 
ceeds in a standard way. Since the computation of the QCD corrections in 
Eq. (|2|) is a one-scale problem (dependent on the ratio Q 2 /m 2 ) the master 
integrals can be expressed in terms of harmonic polylogarithms (HPLs) [7]. 
For the two-loop corrections to 0(e) HPLs up to weight w = 5 are needed, 
H mi ,...,m w {x)i rrij = 0,±1, which depend on the conformal variable x, 

VQ 2 + 4m 2 - y/Q 2 " , Q , 
x = = f= , (3) 



VQ 2 + 4m 2 + y/Ql 

with Q 2 > and < x < 1. The analytic computation of the master 
integrals has used the method of differential equations together with two 
independent numerical checks at selected phase space points based either 
on sector decomposition or the evaluation of a Mellin-Barnes representation. 

The ultraviolet divergences arising in the computation of the Feynman 
diagrams do require renormalization of the strong coupling constant a s , 
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the external (heavy-quark) wave function ip and the heavy-quark mass m 
(typically taken to be the pole mass). The renormalized form factors T\ and 
.T-2 are obtained by adding the appropriate counter-terms CT\ and CTi- The 
renormalization of the bare quantities is performed multiplicatively with all 
necessary constants in the MS-scheme (respectively in the on-shell scheme 
for m) being known [1,8,9]. 

The two-loop corrections to T\ and Ti to 0(e) have been obtained 
in [1] and all terms of 0(e) are new results of [10] . Up to order 0(e°) at 
two-loops full agreement between [1] and the first independent check [10] 
has been established. 

For space-like kinematics the results for T\ and Ti are real whereas for 
time-like kinematics above production threshold (Q 2 > 4m 2 ) both T\ and 
Ti develop an imaginary part, 

^ = 5R^+i9^, (4) 

which can be obtained by means of a suitable (complex) continuation of Q 2 , 
giving rise to x — * —x in Eq. ([2D and all arguments of HPLs. 

3. Exponentiation 

The exponentiation of the heavy-quark form factor is based on the uni- 
versality of soft and collinear radiation and the respective singular terms 
in the electric form factor T\, i.e. the (soft gluon) poles in e and the loga- 
rithms in the heavy-quark mass. The logarithms emerge in the (space-like) 
high-energy limit, Q 2 3> m 2 , i.e. x — ► in Eq. ([5]), as 

L = to (S)- < 5 > 

The magnetic form factor JF 2 is power suppressed by m 2 /Q 2 at high energies. 
The large logarithms of Sudakov type [11] in Eq.([5]) can be resummed. More 
generally T\ fulfills the following evolution equation, see e.g. [4,12], 

- /l2 ^ ln - Fl (^'^' as ' e ) = k K (y >"•><) + k G fy> a '> e ) ' (6) 

where the dependence on the various scales has been made explicit. QCD 
factorization at the scale /i allows to separate to logarithmic accuracy in the 
limit Q 2 ^> m 2 the dependence on the hard scale Q (associated with the 
function G) from that on the heavy-quark mass m (resting in the function 
K). Both functions, G and K, are subject to renormalization group equa- 
tions [12] governed by the same (well-known) cusp anomalous dimension 
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A [13-15], 

In dimensional regularization, the solution of the evolution equation ([6]) pro- 
ceeds analogous to the construction for the massless form factor [6, 16-18]. 
However, in contrast to the massless case where K is a pure infrared counter 
term (dependent on the cusp anomalous dimension A), the heavy-quark 
mass in Eq. ([6]) acts as an additional regulator in the collinear limit. Hence 
the infrared sector differs and K takes particular values in the perturbative 
expansion, see [4,10]. 

The (ultraviolet) renormalized form factor XtiFi reads in space-like kine- 
matics as a function of the renormalized coupling a s (/j, 2 ) 




On the right hand side all quantities are defined in D dimensions as func- 
tions of the Z)-dimensional strong coupling a (to be expanded in terms of 
the ordinary coupling a s in four dimensions). The exponentiation needs an 
additional integration constant C(a(fi 2 , e), e), which is determined through 
matching to fixed-order results in perturbation theory. The perturbative 
coefficients of the function G agree completely with those of the form factor 
for massless quarks. In the latter case the available three-loop informa- 
tion [6, 17, 19] allows to make partial predictions of the high-energy expan- 
sion of T\ based on Eq. (JSj) beyond the two-loop result of [1, 10]. 



4. QCD amplitudes 

For on-shell amplitudes of n-parton processes the factorization ansatz 
has to be extended in comparison to Eq. ©. In D-dimensions a generic n- 
parton amplitude (massless and massive alike) can be expressed as a product 
of functions Jp, <S p and 7i p each governing different regions of kinematics 
(see [20-22]). 

\M P ) = (9) 
Jp ( ^-,a s (/^ 2 ),e j <S P i{ki}, ^2,a s (^ 2 ),ej \H P ) , 
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where the ket-notation |.M P ) implies vectors (and matrices) in color space. 
The jet function J v collects all collinear contributions of the proccess p. It 
is color-diagonal and given as the product of functions J7jji for each external 
parton i. The functions Ju\ are conventionally chosen to be the square 
root of the respective (massless or massive) form factor T\ ■ Soft radiation 
arising from the overall color flow is summarized by the soft function S p , 
which is a matrix in color space and the short-distance dynamics of the hard 
scattering are contained in the hard function TC p , which is to leading order 
proportional to the Born amplitude. 

Eq. ([9]) leads to a remarkably simple relation between n-parton ampli- 
tudes with massless and massive external partons in the small-mass limit. 
If to is small compared to all scales the entire non-trivial mass dependence 
resides in the respective expression for the J™ (i.e. the form factor) so that 
QCD factorization gives [4,5] (see for a review [23]), 



i 





H m) = n (4r l0) ) 3 x H m=o) ' (io) 

i£ {all logs} 

where we the upper index denotes (non-)zero parton masses. The universal 
and manifestly process-independent multiplicative function Z^ ^ is directly 
given in terms of ratios of the respective form factors, 

zt W (^e) = (11) 



j-(m=0) 



where the index i denotes the (massive) parton and a s = a s (/x 2 ). Z^™^ is a 

function of the (process-independent) ratio of scales /j, 2 /m 2 as the (process- 
dependent) scale Q cancels completely in the ratio. This definition however, 
excludes terms with explicit dependence on the number of heavy quarks in 
Eq. (fTTj) (and, therefore in Eq. (fTOl) as well), because at two-loops (and 
beyond) additional process-dependent terms enter in their description [5]. 

Within QCD Eq. (jlOD has been applied in the computation of the (vir- 
tual) two-loop QCD corrections to hadronic top-quark pair-production [24, 
25] in the limit when all Mandelstam invariants are large, s, \t\, \u\ S> to 2 . 
Recent advances beyond the small mass limit have been concerned with the 
exact evaluation of Eq. ([9]) for all kinematics. To that end, in particular the 
soft function 5 p and the structure of its associated anomalous dimensions 
r p have been subject of intense investigation. For the latter quantity T p all- 
order expressions in the case of massless [26, 27] and massive partons 
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have been suggested and explicit results for the two-loop soft anomalous 
dimension of massive partons have appeared recently [29-31]. Based on 
Eq. (jlOp this knowledge has been applied to derive the singularity struc- 
ture of the two-loop virtual amplitude for top-quark hadro-production with 
complete kinematics dependence [32] which was found to be in agreement 
with the small mass limit [24,25]. 



5. Conclusions 

We have provided a brief review of the known perturbative QCD cor- 
rections to the heavy-quark form factor. Thanks to [1] both the electric 
and the magnetic form factor J-\ and J~i are known to two-loop order, an 
important result which has been recently confirmed and extended to higher 
order in e by [10]. The heavy-quark form factor will also be part of future 
calculations of QCD radiative corrections to n-parton scattering processes 
at two-loop order. 

It is a long known fact that the logarithms in the heavy-quark mass m 
in the form factor do exponentiate. These logarithms emerge in the limit 
of small masses upon neglecting all terms m 2 <C Q 2 and within dimensional 
regularization the exponentiation extends to the complete singularity struc- 
ture, i.e. the soft poles in e. The exponentiation is governed by the cusp 
anomalous dimension A as well as the functions K and G. Strikingly, the 
latter is identical for the case of massless and massive quarks, while the 
infrared sector (K) of course differs. 

The knowledge of the complete singularity structure of scattering ampli- 
tudes and, in particular the factorization of amplitudes in soft and collinear 
limits has numerous applications. The factorization property may even be 
used as a tool in practical calculations, a prominent example being the com- 
putation of the virtual corrections to heavy-quark hadro-production at two 
loops in QCD. 
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